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ABSTRACT 

We predict the existence and observational signatures of a new class of objects that assembled 
early, during the first billion years of cosmic time: Obese Black-hole Galaxies (OBGs). OBGs 
are objects in which the mass of the central black hole initially exceeds that of the stellar 
component of the host galaxy, and the luminosity from black-hole accretion dominates the 
starlight. Conventional wisdom dictates that the first galaxies light up with the formation of 
the first stars; we show here that, in fact, there could exist a population of astrophysical objects 
in which this is not the case. From a cosmological simulation, we demonstrate that there are 
sites where star formation is initially inhibited and direct-collapse black holes (DCBHs) form 
due to the photo-dissociating effect of Lyman-Werner radiation on molecular hydrogen. We 
show that the formation of OBGs is inevitable, because the probability of finding the required 
extra-galactic environment and the right physical conditions in a halo conducive to DCBH 
formation is quite high in the early universe. We estimate an OBG number density of 0.009 
Mpc -3 at z ~ 8 and 0.03 Mpc~ 3 at z ~ 6. Extrapolating from our simulation volume, we 
infer that the most luminous quasars detected at z ^ 6 likely transited through an earlier 
OBG phase. Following the growth history of DCBHs and their host galaxies in an evolving 
dark matter halo shows that these primordial galaxies start off with an over-massive BH and 
acquire their stellar component from subsequent merging as well as in-situ star formation. In 
doing so, they inevitably go through an OBG phase dominated by the accretion luminosity at 
the Eddington rate or below, released from the growing BH. The OBG phase is characterised 
by an ultra-violet (UV) spectrum f\ oc A* 3 with slope of /3 ~ —2.3 and the absence of a 
Balmer Break. OBGs should also be spatially unresolved, and are expected to be brighter than 
the majority of known high-redshift galaxies. They could also display broad high-excitation 
emission lines, as expected from Type-I active galactic nuclei (AGN), although the strength 
of lines such as NV and CIV will obviously depend on the chemical enrichment of the host 
galaxy. OBGs could potentially be revealed via Hubble Space Telescope (HST) follow-up 
imaging of samples of brighter Lyman-break galaxies provided by wide-area ground-based 
surveys such as UltraVISTA, and should be easily uncovered and studied with instruments 
aboard the James Webb Space Telescope (JWST). The discovery and characterization of OBGs 
would provide important insights into the formation of the first black-holes, and their influence 
on early galaxy formation. 
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1 INTRODUCTION 

It is now well-established that most present-day galaxies harbour 
a quiescent super-massive black hole (SMBH), with a mass ap- 
proximately one thousand t h of the mass of st ars in the bulge 
i lFerrarese & MerritrJ bOOOt IHaring & Rixl |2004 . Such a corre- 
lation is strongly suggestive of coupled growth of the SMBH 
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and the stellar component, likely via regulation of the gas sup- 
ply i n galactic nuclei from the earli e st times (e.g . ISilk & Reesl 
1 19981 ; lHaehnelt & Kauffmand l200d: iFabianetaU 1 20021: iKinj 
20031 : i Thompson et alJl2005l:lRobertson et alj|200r3: iHopkins et alj 
20091: iNatarai an & Treisterll2009t iTreister et al.ll201ll see however 
Hirschmann et al !|201Cb. 

Since the same gas reservoir fuels star formation and feeds the 
black hole, a connection between these two astrophysical processes 
regulated by the evolving gravitational potential of the dark matter 
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halo is arguably expected. However, understanding when and how 
this interplay commences has been both a theoretical and observa- 
tional challenge for current theories of structure formation. 

In this letter, we explore the formation and evolution of the 
first massive black-hole seeds and the first stars during the earli- 
est epochs in order to explore the onset of coupling between the 
black hole and the stellar component. Our calculation incorporates 
two new physical processes that have only been recently recog- 
nised as critical to understanding the fate of collapsing gas in the 
early universe. The first is the computation of the Lyman-Werner 
(LW) radiation (11.2 — 13.6 eV) that impacts gas collapse in the 
first dark-matter haloes (as it is able to efficiently dissociate the 
H2 molecules, thereby p reventing cooling via molecular hydrogen; 
e.g. lHaiman et al.l d2000l) ). The second is the implementation of our 
growing understanding of the role of the angular momentum of 
the baryonic gas in the co llapse process dLodato & Nataraianl2006l ; 
iDavis & Nataraianl2010l) . Including these two processes within the 
context of the standard paradigm of structure formation, predicts 
the possible existence of a new class of object in the high-redshift 
Universe in which black-hole growth commences before, and con- 
tinues to lead the build up of the galaxy stellar population for a sig- 
nificant period of time. We define an OBG as a phase in a galaxy's 
evolution where post DCBH formation, the BH at least initially 
dominates over the stellar mass and is accreting at a rate sufficiently 
high enough to outshine the stellar component. 

OBGs may provide a natural stage of early black-hole/galaxy 
evolution en route to the most luminous quasars already ob- 
served to be in place at z ~ 6 with estimated black-hole masses 
Mbh — 10 9 Mq. Observationally, such objects should appear sim- 
ilar to moderate-luminosity AGN, but with very low-luminosity 
host galaxies, and low metallicities. 



2 METHODOLOGY 

Our model is a modified version of I Agarwal et al.ld2012l) . A12 here 
after, where they identify the sites of DCBH formation by calcu- 
lating the total amount (spatial and mean) of LW radiation seen 
by any given halo within a cosmological N-body dark-matter only 
simulation using a semi-analytic model for the star formation in 
these haloes. The key features of A12 are summarised below: 

(i) The DM only N-body simulation is run from z = 30 to z = 
6 with a box size of ~3.4Mpc h~ . 

(ii) The minimum halo mass resolved with 20 particles is ~ 
10 Mq, an order-of-magnitude lower t han the minimum ha lo 
mass that can host a Pop III star at 2 = 30 dTegmark etail 19971) . 

(iii) Both Pop III and Pop II star formation are allowed, and halo 
histories are tracked in order to determine if a halo is metal free. 

(iv) Many realisations of the model are run to study the effect 
of different LW escape fractions, Pop II star formation efficiencies, 
gas outflow rates due to supernova feedback, number of Pop III 
stars forming per halo and reionisation feedback. The results pre- 
sented here are based on the run with a LW escape fraction of 1.0 
and a Pop II star-formation efficiency of 0.005 (see A12 for more 
details). 

(v) The LW specific intensity, JlwEI is computed self- 
consistently and with two components: a global and a local con- 
tribution. A pristine halo is considered for Pop III star formation or 

1 J combined with a super-script or sub-script expresses the specific inten- 
sity of LW radiation in units of erg s — 1 cm — 2 Hz — 1 sr _1 . 
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Figure 1. The temperature/spin distribution of pristine massive dark-matter 
haloes that are exposed to Jlw ^ >^crit- The virial temperature, T v ; r , is 
plotted against the halo spin parameter, A, for all DCBH candidates with 
Toomre stability parameter Q c = 1.5, 2, 3. The nearly-vertical solid curve 
represents the value of A max and all haloes with spin less than this critical 
value lie to the left of this line (red points), whereas those with larger spin lie 
to the right (black points). The upper limit on the scale-length of the hosted 
disc given the allowed T v ; r — A max combination is marked as the dashed 
line. Any halo in the yellow region, i.e. below the dashed line and to the 
left of the solid vertical curve will host a DCBH (blue points) in our model. 
Note that the yellow region shrinks as Q c decreases, thereby reducing the 
probability of finding a halo that can form a DCBH. Inset: A zoom-in on 
the T v j r — A distribution of the four DCBH candidates in our fiducial case 
with Q c = 3. 

treated as a DCBH candidate depending on the halo's virial tem- 
perature and the Jlw that it is exposed to. 

(vi) The number density of DCBH sites can be up to 0.1 Mpc~ 3 
at z = 6, much higher than previously anticipated foiikstra et al.l 
120081) . 

In the present study we refine the model of DCBH formation 
and also follow the subsequent growth of these seeds as identified in 
A12. We discuss new additions to the A12 model in the subsections 
below. 



2.1 DCBH forming haloes 

A DCBH for ms in our model if a pristine massive halo is exposed to 
J ^ JrrH PI jWolcott-Green et alJ[2oT 11) and satisfies both the spin 
and size criterion required for th e disc to withstand fragmentation 
dLodato & Nataraia nl2006ll2007l . LN06 and LN07 hereafter). Note 
that the J cr i t is always produced by stellar sources external to the 

2 Note that J cr i t is the critical level of LW radiation required by a pristine 
atomic cooling halo to undergo direct collap se. For Pop III stars, J cr jt ~ 
1000 and for Pop II stars J crit ~ 30 - 100 iWolcott-Green et al.l20lH) 
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Figure 2. Predicted redshift evolution of Mbh and M„ for the four OBGs 
in our simulation down to z = 6. We show the tracks for different frac- 
tions of the Eddington accretion rate /i; m = 1.5, 0.75 (fiducial), 0.1. The 
parameter z g denotes the redshift when the galaxy becomes an OBG and 
first appears on this plot. The shaded portion represents the region where 
the BH's accretion rate would have to be larger than the Eddington limit 
for the galaxy to qualify as an OBG. In solid-black l ines, we show the lo cal 
Mbh~ Mbuige relation and the I— a error in the fit |Haring & Rix 20041) . 



pristine DCBH candidate halo as internal sources would pollute the 
gas inside the halo. We assume the collapsing gas in pristine haloes 
will settle into a disc whose stability against fragmentation deter- 
mines whether it will be able to collapse to a BH or will fragment 
into star-forming clumps (LN06, LN07). Assuming that initially the 
baryons have the same specific angular momentum as the halo, the 
halo must have a spin, A, lower than a characteristic value, A max , 
for a given Toomre stability parameter Q c , for which the pristine 
gaseous disc is exactly marginally stable and above which no ac- 
cretion can take place onto the central region (LN06). The critical 
spin is given as 



Amax — 



«;'d 



(l) 



where md is the disc mass expre ssed as a fixed fr action (0.05) of the 
total baryonic mass in the halo dMo et alii 1998b . jd is the specific 
angular momentum of the disc that is also a fixed fraction (0.05) 
of the halo's overall angular momentum (LN06), T v i r is the virial 
temperature of the halo and T gas is the temperature of the gas in the 
disc which depends on whether atomic or molecular hydrogen is 
the dominant cooling species. In our case, since the halo is exposed 
to J ^ Jcrit, the dominant coolant is atomic hydrogen and T sas is 
set to 8000 K accordingly. The absence of the H2 cooling channel, 
which can cool the gas down to 200 K, is ensured by J cl i t . 

The second condition comes from the limit that, for a given 
channel of gas cooling in the pristine disc (i.e. H or H2), the disc 
must be cooler than a characteristic temperature (proxy for size of 
the disc), expressed as T max , above which the gravitational torques 
required to redistribute the angular momentum become too large 
and can disrupt the disc. 



T max -T gas ^— 1 + A/Bn/mdM 



2/3 



(2) 



where a c is a dimensionless parameter (0.06) relating the critical 
viscosity to the gravitational torques. This provides a mass estimate 
for the assembling DCBH (LN07): 



8Aj d fT ( 



mlQ c 



gas 

T ■ 

J- vir 
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for A < A max . 

We find that the inclusion of these two criteria for efficient 
angular momentum transport and accretion within the disc, in ad- 
dition to our existing framework for the treatment of LW radiation 
feedback, fundamentally alters the progression of structure forma- 
tion in these haloes, and impacts the observable characteristics of 
stars and the BHs within them. 

We plot the T v i r — A distribution of pristine atomic cool- 
ing haloes that are exposed to Jlw ^ Jcrit, for different values 
of Q c in Fig. 1. The haloes with spin in the range A < A max 
are marked in red, with the almost-vertical solid curve represent- 
ing A max , whereas the ones with A > A max are plotted in black. 
The size constraint, in order for the disc to withstand fragmenta- 
tion, is denoted by the dashed curve. These limits together, marked 
by the yellow region, constrain DCBH formation to a small al- 
lowed domain in the T v i r — A plane. This still yields copious num- 
bers of DCBHs (blue points) with a co-moving number density of 
0.03 Mpc -3 in our fiducial case with Q c — 3, and fu m = 0.75 
(see the following section for /n m ). Note that LN06, LN07 require 
the gas disc to be marginally stable/unstable i.e. Q c ~ 1 — 2. Given 
that the actual disc parameters are uncertain, we choose Q c — 3, 
which allows us to include most of the DCBH candidates for rea- 
sonable disc parameters, and for which the disc sizes are not too 
large. 

2.2 Growth of a DCBH 

Haloes hosting DCBHs are initially not massive enough and are 
not polluted enough to lead to Pop II star formation (Schneider et 
al. 2002). In our model, once a DCBH forms, the halo may host 
metal-free Pop III stars, as it is plausible that the gas disc could un- 
dergo fragmentation post DCBH formation. Following this epoch, 
however, the BH grows by accreting gas available in the halo, un- 
challenged by any further star formation. Note that Pop III stars 
forming in a halo could in principle pollute it if the halo is mas- 
sive enough to retain the metals. We thereby set a mass threshold 
of M > 10 s Mq for polluted haloes to form Pop II stars. It is rea- 
sonable to assume that prior to this, the gas reservoir is still massive 
enough to feed the central BH. At this stage the accreting DCBH 
might appear as a mini-quasar, but of essentially zero metallicity. 
Once the halo is massive enough (A/dm > 1O 8 M0), due to the 
metal pollution by the first generation of Pop III stars, the gas can 
ultimately cool leading to the formation of Pop II stars. The stellar 
component and the BH from this point on begin to grow in tandem, 
marking the onset of the OBG phase. How the two components 
evolve in detail is sensitive to the accretion rate and subsequent 
merging history - a parameter space that we have explored exten- 
sively. 

We assume in our model that a DM halo is initially comprised 
of hot gas, Mhot = /bA/DM, where /b is the universal baryon frac- 
tion and A/dm is the halo's current DM mass. Mhot transitions to 
the cold gas phase, M co id, by collapsing over the dynamical time of 
the halo, tdyrEl Further details on Pop II star formation and super- 
nova feedback are given in A12. We do not treat metal diffusion 
explicitly from external stellar sources. This could, in principle, 



tdyn is defined as the ratio of the halo's virial radius to the circular ve- 
locity defined for the infall mass and infall redshift 
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pollute a pristine halo. However, the metal diffusion through the 
inter-galactic mediu m is highly inefficient ll Johnson et al.ll20~12al : 
iMuratov et al]|2012h . and even if metals do reach a pristine halo, 
they need to diffuse within and make their way to the central region, 
and t his process has a larger typical timescale dCen & Riquelmej 
l2008h than DCBH formation. 

Once a DCBH forms in a halo, it is allowed to grow at a fixed 
fraction /n m of the Eddington accretion rate, with the radiative effi- 
ciency set to 10 %, assuming that both the cold gas and hot gas can 
be accreted by the BH. The upper limit of the accretion rate is set 
by the parameter /ii, n that we vary between individual runs. If the 
total gas available during our integration time steps for accretion 
is less than this fraction, the total mass available sets the accretion 
rate. We run our model for /n m = [1.5, 1.0, 0.75, 0.1] to explore 
the parameter space. Our fiducial case corresponds to a LW escape 
fraction of 1.0, Pop II star formation efficiency of 0.005, Q c — 3 
and /n m = 0.75. In haloes which host a DCBH but are not massive 
enough to form Pop II stars, the gas is assumed to be hot and dif- 
fuse (i.e. has not condensed over the dynamical time of the halo). In 
haloes which host a DCBH and a Pop II stellar component, both the 
hot and cold phases of gas are assumed to contribute to the accre- 
tion process. The total mass accreted by the DCBH is split into hot 
and cold components depending on the ratio of the hot and the cold 
gas reservoirs. We do not assume any feedback from the accreting 
DCBH affecting star formation in the galaxy. We do this to avoid 
inserting a correlation between the BH and stars by assuming such 
a feedback loop since the precise nature of accretion and feedback 
in galactic nuclei is largely unknown at z > 6. 

To summarise, the total gas mass available for accretion, and 
hence the total mass accreted by the DCBH, depends on whether 
Pop II stars are forming in the halo or not. If there is no assembling 
Pop II stellar component, the DCBH is assumed to accrete from the 
hot gas reservoir. If the halo hosts a Pop II stellar component, the 
amount of gas accreted is split according to the ratio of available 
hot and cold gas. 



3 RESULTS AND DISCUSSION 

For Q c = 3, we find four OBGs, named 01-04, in our simula- 
tion box. The stellar and black-hole growth tracks of these OBGs 
are shown in Fig. 2, colour coded as Ol-blue, 02-red, 03-green, 
04- purple, respectively. The fiducial case (/n m = 0.75) is marked 
by the solid curves and the open triangles and circles denote the 
time-steps at fn m = 0.1, 1.5. The grey shaded region is where the 
BH would need to accrete at super-Eddington rates for the galaxy 
to appear as an OBG. Since these objects have Mbh > A/*, the 
nuclear emission can dominate the starlight even when accreting at 
significantly sub-Eddington rates in the non-shaded region. 

3.1 Observational predictions 

Inclusion of these physical processes enables the calculation of 
the observational signatures of OBGs. Haloes that harbour grow- 
ing DCBH seeds with no (or little) associated stellar component 
merge into haloes that have formed the first and second genera- 
tion of stars. We compute the observed spectral energy distribution 
(SED) of these copiously accreting DCBH seeds and the population 
of stars within the OBGs. To model the stellar component of the 
SED, we use the stellar masses and ages fro m our merger tree an d 
derive the spectrum using STARBURST99 dLeitherer et all 1 19991) . 
The accretion disc spectrum is modelled as a radiating blackbody 
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Figure 3. OBG candidates and their observability with JWST at z ~ 6. We 
plot the observed flux density in AB magnitudes for all the OBGs, while 
varying the maximal allowed accretion rate, /i; m . When the BH accretion 
dominates the spectrum there is virtually no Balmer break, and the UV slope 
is obviously fixed by the accreting black hole. The black points denote the 
flux limits and bandpass widths of NIRCam (circles) and MIRI (squares) 
wide filters, assuming a 10,000 second exposure with JWST and a S/N ratio 
of 10. The shaded area marks the wavelength region shortward Ly«, where 
intergalactic neutral hydrogen is expected to completely absorb the OBG 
signal. 



with a tempera ture profile of the disc given by the alpha-disc model 
jPringld[T98ll) . Since OBGs are expected to lie in haloes of low 
metallicity we do not include any dust absorption in our models. 
Note that an OBG is characterised by possessing an actively accret- 
ing BH and an underlying stellar population that is Pop III or Pop 
II or both, however an OBG might appear as what has often been 
referred to as mini-quasars in the very early stages of its evolution 
when no stellar component is found in the DCBH host galaxy. 

The UV-optical SED of an OBG is inevitably dominated by 
the accretion onto the central black hole. The predicted SEDs of 
the OBGs over the wavelength range observable with the NIRCAM 
(Near Infra-red Camera) and MIRI (Mid-Infra red Instrument) in- 
struments aboard NASA's proposed James Webb Space Telescope 
(JWST) are shown in Fig. [3] The stellar spectrum dominates the 
OBG spectrum only when the BH is limited to < 0.1 f c( jd at all 
times. However, such a low rate would be incompatible with the 
dramatic mass growth rates expect ed of the most massive, early 
black holes (e.g lsiiacki et alj2009l) . 

We note that the magnitude of our brightest OBGs could be 
mAB ~ 25, comparable to the brightest putative Lyman-break 
galaxies uncovered at z ~ 7 in ground-based surveys such as Ultra- 
VISTA l lBowler et alj2012h . However, it will be hard to distinguish 
OBGs from Lyman-break galaxies with ground-based imaging be- 
cause, while the predicted UV continuum black-hole emission is 
expected to be relatively blue (with a UV slope j3 ~ —2.3, where 
/a <x A' 3 ), it is not significantly bluer than that dis played by the 
gener al galaxy population at these early times (e.g. iDunlop et alj 
|2012|) . OBGs are also of course expected to display a negligible 
Balmer break. However, while it is interesting that the stack of 
the brightest Lyman-break galaxies in Ultra VISTA shows at most a 
very weak Balmer break, it will still require extremely high-quality 
mid-infrared photometry to conclusively rule out the possibility 
that the UV-optical SED of a putative Lyman-break galaxy is in- 
compatible with that produced by a very young stellar population. 
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Figure 4. Size versus magnitude relation. The HUDF galaxies, not cor- 
rected for point spread function (PSF), at z = 6, 7, 8 are represented by up- 
right triangles, downward triangles and squares respectively. Observational 
limits and the PSF of the NIRCam are plotted as the straight lines. The 
OBGs 01-04 in our sample, denoted by the arrows pointing left, would 
be unresolved objects that could be brighter than the galaxies. Note that we 
have excluded the fn m = 0.1 case from this plot as the mAB for the OBGs 
is quite high. 



Identification of OBGs amidst the observed high-redshift 
galaxy population will therefore require high-resolution imaging 
with HST and ultimately JWST. With A/ BH > M„ , an OBG accret- 
ing at a reasonable fraction of the Eddington limit should certainly 
appear unresolved and point-like in high-resolution rest- frame UV 
(observed near-infrared) imaging. As illustrated in Fig. 4, it is al- 
ready known that the vast majority of faint high-reds hift galaxies 
uncovered via deep H ST imaging are resolved (see lOesch et all 
l20ld : IOnoet~ai]|2012l) . but this does not rule out the existence of 
an OBG population with a surface density < 1 arcmin~ 2 , and HST 
follow-up of the brighter and rarer high-redshift (z — 7) objects 
uncovered by the near-infrared ground-based surveys is required to 
established whether or not they are dominated by central black-hole 
emission. 

Finally, like all active galactic nuclei, OBGs are expected 
to display broad-line emission from highly excited species in the 
vicinity of the black hole. However, as OBGs have very low metal- 
licity, it is unclear whether lines such as NV and CIV are expected 
to be detectable even given high-quality near-infrared spectroscopy, 
and Lyman-a is often severely quenched by neutral Hydrogen as 
we enter the epoch of reionisation. 

Thus, the best observational route to establishing whether 
OBGs exist, and if so constraining their number density (and ul- 
timately their evolving luminosity function), appears to be via deep 
imaging of putative high-redshift Lyman-break galaxies, with suf- 
ficient angular reolution to prove they are unresolved, coupled with 
sufficiently accurate photometry to prove any poi nt-like object s 
cannot be dwarf star contaminants (see, for example |Punlop|2012h . 

The discovery of an OBG could in principle settle the long 
standing debate on whe ther DCBHs can form and be the seeds of 
the first SMBHs (A12 iDiikstra et alj[200j; IVolonteri et alj|2008l; 



ACKNOWLEDGEMENTS 

The authors would like to thank Rychard Bouwens for providing 
the HUDF data used in Plot. [4] JSD acknowledges the support of 
the European Research Council via an Advanced Grant, and the 
support of the Royal Society through a Wolfson Research Merit 
award. 



REFERENCES 

Agarwal B., Khochfar S., Johnson J. L., Neistein E., Vecchia 
C. D., Livio M., 2012, MNRAS, 425, 2854 

Bellovary J., Volonteri M., Governato F, Shen S., Quinn T., Wad- 
sleyj., 2011, ApJ, 742, 13 

Bowler R. A. A. et al., 2012, MNRAS, 426, 2772 

Cen R., Riquelme M. A., 2008, ApJ, 674, 644 

Davis A. J., Natarajan P., 2010, MNRAS, 407, 691 

Dijkstra M., Haiman Z., Mesinger A., Wyithe J. S. B., 2008, MN- 
RAS, 391, 1961 

Dunlop J. S., 2012, eprint arXiv,1205,1543 

Dunlop J. S. et al., 2012, eprint arXiv, 1212, 860 

Fabian A. C, Wilman R. J., Crawford C. S., 2002, MNRAS, 329, 
L18 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Haehnelt M. G„ Kauffmann G., 2000, MNRAS, 318, L35 

Haiman Z., Abel T., Rees M. J., 2000, ApJ, 534, 11 

Haring N., Rix H.-W., 2004, ApJ, 604, L89 

Hirschmann M., Khochfar S., Burkert A., Naab T., Genel S., 

Somerville R. S., 2010, MNRAS, 407, 1016 
Hopkins P. F, Murray N., Thompson T. A., 2009, MNRAS, 398, 

303 

Johnson J. L., Vecchia C. D., Khochfar S., 2012a, eprint arXiv, 
1206, 5824 

Johnson J. L., Whalen D. J., Holz D. E., 2012b, eprint arXiv, 1211, 
548 

King A., 2003, ApJ, 596, L27 

Leitherer C. et al., 1999, ApJ, 123, 3 

Lodato G, Natarajan P., 2006, MNRAS, 371, 1813 

Lodato G, Natarajan P., 2007, MNRAS: Letters, 377, L64 

Mo H. J., Mao S., White S. D. M., 1998, MNRAS, 295, 319 

Muratov A. L., Gnedin O. Y., Gnedin N. Y., Zemp M., 2012, eprint 

arXiv, 1212, 909 
Natarajan P., Treister E., 2009, MNRAS, 393, 838 
Oesch P. A. et al., 2010, ApJL, 709, L21 
Ono Y. et al., 2012, arXiv, 1212, 3869 
Pringle J. E., 1981, Ann. Rev. Astron. Astrophys., 19, 137 
Robertson B., Hernquist L., Cox T. J., Matteo T. D., Hopkins P. F, 

Martini P., Springel V, 2006, ApJ, 641, 90 
Sijacki D., Springel V, Haehnelt M. G, 2009, MNRAS, 400, 100 
Silk J., Rees M. J., 1998, A&A, 331, LI 

Tegmark M., Silk J., Rees M. J., Blanchard A., Abel T., Palla F, 

1997, ApJ, 474, 1 
Thompson T. A., Quataert E., Murray N., 2005, ApJ, 630, 167 
Treister E., Schawinski K., Volonteri M., Natarajan P., Gawiser 

E., 2011, Nature, 474, 356 
Volonteri M., Lodato G., Natarajan P., 2008, MNRAS, 383, 1079 
Volonteri M., Natarajan P., 2009, MNRAS, 400, 191 1 
Wolcott-Green J., Haiman Z., Bryan G. L., 2011, MNRAS, 1673 



Volonteri & Nataraianll2009t|johnson et al.ll2012bl ; lBellovary et au 
201 lh . Uncovering this population holds great promise for under- 



standing the onset of black hole and host galaxy growth. 



© 0000 RAS, MNRAS 000,[TJ|5] 



